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Abstract 

We examine a superpotential for an SO(10) GUT and show that if the parameters 
of the superpotential are in a certain region, the SO(10) GUT has an intermediate 
symmetry SU(2) L (g) SU(2) iJ (g) SU(3) C ® U(l) B _ i which breaks down to the group 
of the Standard Model at an intermediate scale 10 10-12 GeV. In the model by the 
breakdown of the symmetry right-handed neutrinos acquire mass of the intermediate 
scale through a renormalizable Yukawa coupling. 
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1 Introduction 



When we construct a Grand Unified Theory(GUT) based on SO(10) in general, we 
have singlet fermions under the Standard Model(SM) -what we call right-handed neutrino. 
Under the SM right-handed neutrinos can have Majorana masses because they are singlet. 
Then the scale of the right-handed neutrinos(= M VR ) is expected to be a scale below which 
the SM is realized. 

It is well known that in the Minimal Supersymmetric Standard Model (MSSM) the 
present experimental values of gauge couplings are successfully unified at a unification 
scale Mu ~ 10 16 GeV 0. This fact implies that if we would like to consider the gauge 
unification, it is favorable that the symmetry of the GUT breaks down to that of the 
SM at the unification scale. In this case the scale of the right-handed neutrinos M VR is 
expected to be the unification scale Mu- This means also there is no intermediate scale 
between the Supersymmetry(SUSY) breaking scale and the unification scale. 

On the other hand it is said that M VR ~ lO lo ~ 12 GeV0. The experimental data on the 
deficit of the solar neutrino can be explained by the Mikheyev-Smirnov-Wolfenstein(MSW) 
solution Q. According to one of the MSW solutions, the mass of the muon neutrino seems 
to be m v 1CT 3 eV. Such a small mass can be led by the seesaw mechanism ||: A muon 
neutrino can acquire a mass of O(10" 3 ) eV if the Majorana mass of the right-handed 
muon neutrino is about 10 12 GeV. 

Then how can the right-handed neutrinos acquire mass of about 10 12 GeV? It was our 
question in our previous paper |J, because if we take the prediction of the MSSM serious, 
M VR is expected to be Mjj ~ 10 16 GeV. Our point of view was that it is more natural 
to consider that one energy scale corresponds to a dynamical phenomenon, for instance a 
symmetry breaking. Mass is given by a renormalizable coupling is also the crucial point 
of our view. This idea is consistent with the survival hypothesis. Thus we were led to a 
possibility that a certain group breaks down to the SM group at the intermediate scale 
at which right-handed neutrinos gain mass through a renormalizable coupling. 

In the previous work we have searched possibilities to construct such a SUSY SO(10) 
GUT with an intermediate symmetry. We have a possibility to construct a SUSY SO(10) 
GUT with an intermediate symmetry SU(2) £ ® SV(2) R <g> U(1) B _ £ <8> SU(3) C (= G2231) Q 
which breaks down to the SM group at an intermediate scale M VR ~ 10 10_12 GeV where a 
right-handed neutrino gains mass. 

In such a scenario, as we showed in the previous work, to make the model consistent 
with the gauge unification we have to introduce several multiplets at the intermediate 
region between the GUT scale and the intermediate scale, in addition to ordinary matters, 
three generations of quarks and leptons and a pair of so-called Higgs doublets. 

Although we showed a possibility to construct a SUSY SO(10) GUT with an interme- 
diate symmetry G2231 it is not trivial whether it is actually possible to construct such a 
GUT since there are many extra fields in the intermediate region. We did not show the 
superpotential for the theory explicitly which can realize such a scenario that we have 
suggested in ref.0. 

1 We use a notation G/ m „... to represent SU(1) ® SU(m) ® SU(n)....If I = 1, it means U(l). 
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The purpose of this paper is to show an explicit form of a superpotential for a SUSY 
SO(10) GUT to construct a SUSY SO(10) GUT whose symmetry breaks down to G2231 
at a GUT scale Mjj and G2231 breaks down to the SM symmetry at the intermediate scale 
M VR . 

We give the scenario and the model briefly in sect. 2 where we give a candidate for 
the matter content in the intermediate region (the spectrum (]!])). Then in sect. 3 we 
show the most general form of the superpotential and a symmetry breaking condition 
as preparation for our analysis. In sect. 4 first we calculate parameters of the theory, 
namely parameters appearing in the superpotential, which produce the spectrum ([]]) at 
the intermediate region. Then we show the exact parameters which realize the MSSM 
below M VR . Finally (in sect. 5) we give a summary and discussion. 



2 Scenario and Model 
2.1 Scenario 

We construct a SUSY SO(10) GUT whose symmetry breaks down to G2231 at a GUT scale 
M\j and G2231 breaks down to the SM symmetry at the intermediate scale M VR . When 
G2231 breaks down to the SM symmetry the right-handed neutrinos gain mass through a 
renormalizable Yukawa coupling. 

Let us first recapitulate the content of the previous work0. To achieve the gauge 
unification in the scenario we have to introduce a certain combination of multiplets. 
Because in our model right-handed neutrinos acquire mass of 0(M UR ) via a renormalizable 
Yukawa coupling by the symmetry breaking, we have to introduce at least a pair of 
(1,3,1,6) + h.c multiplet under (j2231- We adopt the normalization for U(1)b~l, T} b = 
diag(— 1, —1, —1,3). When we introduce only (1,3,1,6) + h.c multiplet in addition to the 
ordinary matter, gauge couplings do not unify. Then we have to introduce certain matter 
content under £2231- 

We found very many candidates for matter content in the intermediate region between 
the GUT scale and the intermediate scale which lead the gauge unification. Among them 
we showed two candidates for the matter content as the simplest example. In this article 
we use another candidate which was not showed in the previous paper. In the examples 
appearing in it a (1,3,1,0) multiplet under G2231 was not included. In constructing a GUT 
following the idea, however, we have to introduce a (1,3,1,0) multiplet in the intermediate 
region. The reason why we have to introduce a (1,3,1,0) multiplet is stated in the appendix 
A. Thus we have to use another candidate for matter content. 

The matter content other than quarks and leptons (including right-handed neutrinos), 
which we assume survive until G2231 breaks down to the SM group at the intermediate 
scale, are given below. 
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(1,3,1,-6) 1 (1,3,1,6) 1 responsible fori/Rinass 
(2, 2,1,0) 2 ordinary Higgs doublets 

(2,1,3,-1) 1 (2,1,3,1) 1 

(2,1,1,3) 1 (2,1,1,-3) 1 1 ] 

(1,3,1,0) 1 
(1,1,8,0) 1 

In this list, for example, (1,3,1,-6) 1 stands for that the representation of the matter 
under G2231 is (1,3,1,-6) and its number is one. When we have the particle content listed 
here in the intermediate region the unified coupling au(Mjj) is about 1/18 if we take the 
intermediate scale to be 10 12 GeV. As a candidate which contains (1,3,1,0), this candidate 
leads the smallest unified coupling. 

In our scenario, at the GUT scale Mjj where SO (10) breaks down to G2231 almost of 
all particles have mass of O(Mu) while the particles listed in (|I|) as well as quarks and 
leptons are left massless. Then at the intermediate scale where G2231 breaks down to G231 
the SM group all extra multiplets in (|1|) besides a pair of Higgs doublets and right-handed 
neutrinos have mass of 0(M UR ), that is, they decouple from the spectrum. Thus below 
M VR the MSSM is realized. 

2.2 Model 

2.2.1 Matter content 

To have multiplets (H) and quarks/leptons at the intermediate region we introduce fol- 
lowing multiplets of SO(10). 





SO(10) 


C2231 


H 


10 


(2,2,1,0),... 


A 


45 


(1,3, 1,0), (1,1, 8,0),. 


$ 


126 


(1,3, 1,-6), (2, 2, 1,0) 


$ 


126 


(1,3, 1,6), (2, 2, 1,0),. 


A 


210 


(1,3, 1,0), (1,1,8,0),. 


i=l~4 


16 


(2, 1,3,-1), (2, 1,1, 3) 




16 


(2, 1,3,1), (2, 1,1, -3) 



(2) 



quarks/leptons 



In this list numbers in the column of SO(10) means SO(10) representation. In the last 
column we show what representation in (H) is contained in the corresponding SO (10) 
multiplet. 

By the requirement that the right-handed neutrinos get mass through a renormalizable 
coupling, we introduce 126 and 126. As a candidate of ordinary Higgs doublets 10 is 
introduced. There are other candidates for ordinary Higgs doublets in 126 and 126. Then 
the ordinary Higgs doublets will be a mixture of these three. To break SO(10) to the SM 
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group via G2231, namely to have the intermediate symmetry G2231, we use 45 and 2100. 
These also contain (1,3,1,0) and (1,1,8,0). 4 16's and 1 16 represent 4 generation + 1 
anti-generation. The reason why we introduce a pair of 16 and 16 is that they contain 
(2,1,3,-1) + h.c and (2,1,1,3) + h.c . 

At this stage the matter content (||D is just a candidate which may realize our scenario. 

As we will see, we can write down the superpotential with these matter which realize 
our idea. 

2.2.2 Singlets under the SM group 

In the SO(10) multiplets @ there are many singlets under the SM symmetry (see appendix 
B for the meaning of subscripts 1,...,0): 



Field : Component Little Group 

A : &12+34+56 = a. G2231 

: a 78+90 = P G24I 

'■ 01-2i,3-4i,5-6i,7-8t,9-Oi = SU(5) 

^ : 01+2i,3+4i,5+6i,7+8i,9+Oi = SU(5) 

A : ^7890 = a G224 

^1234+3456+5612 = b G223I 

^(12+34+56)(78+90) = c G23H 

^=1^4 : VW~4 SU(5) 

* : if) SU(5) 



where a,b,... stand for vacuum expectation values (VEV) of the corresponding fields. 
Little group means a remaining symmetry when only a corresponding component has a 
VEV. For example, when only a gets a VEV SO(10) breaks down to GW- 

Among them a, b and a are G2231 singlets and hence their order of magnitudes is 
expected to be the GUT scale My ~ 10 16 GeV. By assumption that SO(10) breaks down 
to G2231 at the GUT scale, b or a must be of order M v . Others must be of order at most 
M UR = Mjje by assumption because they are not G2231 singlets. Also <p is required to be 
of order M„ K , 

4~M VB (=Mue) (4) 

because it gives masses to the right-handed neutrinos. Of course, as we will see, there 
are constraints among VEVs in addition to the well known constraints - F-flat and D-flat 
condition because we require certain multiplets must have mass of 0(M„ fl ). 

2 Using only 210 it is impossible to break SO(10) to G231 through G2231 0- We can break SO(10) to 
the SM group via G2231 using 45 + 54. In this case if there is no multiplet which contains (1,3,1,0) other 
than 45 (3,1,1,0) is also massless. The reason is that mass terms for (1,3,1,0) and (3,1,1,0) come from 
the mass term of 45 and the vacuum expectation value of 54 through the coupling 45 2 54 and because of 
D parity^] they are same as each other's. Thus we can get rid of the possibility of using 45+54. 
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3 Preparation 



3.1 Superpotential 

With the multiplets (0) the most general form of the superpotential W is written as 

W = W mass + W mt + W*. (5) 
W mass consists of the most general bilinear terms: 

W mass = ^M H H 2 + M $ $$ + ^M A A 2 + l -M A A 2 + M^ 4 . (6) 

We define only \l/ 4 has a mass term with because by a redefinition of ^4 , namely 
by a rotation among \J/ i=1 _ 4 , it is possible that only the new ^4 has a mass term with 
We require all mass parameters are 0{Mu) because Mjj is the natural order for them. 
Wi n t has the most general interaction terms without 16 and 16: 



W mt = Y H * A H$A + Y Hl5A H$A + -Y A A 3 + Y $A $A$ + Y^A® 



1 

3! 



+ -Y AA 2A 2 A + -Y A 2 A AA 2 . (7) 

We require all Yukawa couplings are at most 0(1). More exactly, as an expansion pa- 
rameter for the perturbation we require they are at most 0(1). As a expansion parameter 
for the perturbation they appear with multiplied by a certain overall factor. The overall 
factors for each couplings are given in appendix B.3. 

Finally, represents the most general interaction terms with 16 and 16. 

4 4 

1=3 i=2_ ij (8) 

+ '£ i y ij * i * j H + y , **H. 

ij 

By the same reason that only \l/ 4 has a mass term with only ^3 4 have couplings 
with A and only ^2,3,4 have couplings with A. 

To see in which direction the gauge group SO(10) can break down we examine the 
D-term and the F-term conditions. 

3.2 D-flat condition 

To keep the SUSY all D-terms must be zero up to SUSY braking scale: 

$t T |$ + + J- q\T^i + ^Tpr + A f T£A + A ] T%A = 0. 



Since the D-term for real representations automatically vanishes [pi 10 1, 
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2(I0I 2 -I0l 2 ) + (EI^I 2 -H 2 ) = o (9) 

1=1 

must be satisfied. The factor 2 reflects the difference of U(l) charge which corresponds 
to a broken generator. 

Later we put ^'s and ip zero. In this case 

|0| 2 -|0| 2 = O. (10) 

3.3 F-flat condition 

First we examine the F-flat condition for 16 and 16 with a mass term for (1, 2, 1, —3) + h.c 
component because the singlet components of 16 and 16 are contained in it and therefore 
there is a relation between the mass term and the F-flat condition. By such an examination 
we see that ipi and ip should be zero though it is not a strict reason for it. 

The F-flat condition for 16 and 16 are as follows: (See appendix B to know how to 
calculate Clebsch-Gordan (CG) coefficient) 

dW - 

— = 2^^,0 = 0, (11) 
dW 4 — — 

— = 2 £ y 2j ij j( p - Y 9A2 (V6ia + 2i/?)V = 0, (12) 

— — = 2 V3jipj4> - Y 9A3 (\/6ia + 2ip)ip 

- y*A3(2"\/6a + 6v^6+12c)^ (13) 
= 0, 



dW 



4 

j'=i 

- y^ A4 (2>/6a + 6v^6 + 12c)^ + M^ (14) 
= 0, 



dW — 4 

-==- = 2 J /V0 + E- y »^( v ^ a + 2 W< 

i=2 

4 

• ^n,Ai(2v / 6a + 6v / 2fe + 12c)^ + M^ 4 (15) 

: 0. 



dip 
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By the way in the intermediate region where G2231 is realized, (3 = c = and the mass 
term for (l,2,l,-3)+h.c is given by 

/ \ 



d 2 W 

dipidtp 



— \f6iYipA20t 
-\/6iYy A3 a - 2v / 65 / * A3 (a + \/3b) 



(16) 



If (j),(j),ipi,ip ~ 0(e), using F-flat conditions (0) - (|14D, all of elements of the mass term 
for (l,2,l,-3)+h.c, fllBD , are calculated to be of order M VR . This, however, contradicts with 
the mass spectrum (p]). Though we may be able to make the some elements of the mass 
term 0(M[/), for example, by making ip 0(e 2 ) (with an appropriate value of ipi,<fi ~ 0(e)), 
we put ipi and if> zero since what we try to do is to show a possibility of SUSY SO (10) 
GUT with an intermediate scale and to take if>i = if) = as the solution of the F-flat 
conditions for 16 and 16 is the easiest way of it. 
Then other F-term conditions follows: 



dW 



oa Fo "V h Yaa2(3 1 YaC 1 u ■ Y * A ^ n 



;i7) 



dW 
~db 



dW 
dc 



OA /O'V ■ Fa&2 

= 24 v 2 z ia 2 a aa f= 1 7= 

3V2 18 y/2 

+ 24v / 2^1a2a/5c+^= + M a 6+ y "' A 



18 



10 V2 



0. 



16 v 6 2 Ya 2 a a c + — — j=- + Mac + 



9^ 



10 



0. 



(19) 



dW 

da 



oa fo ' v h ^ Y AA^ab Y AA 2 f3 c 
= 24 V2 1 Y A 2 A ab -^=— 

o fa -V 2 . w , V^^a^^ 

+ 8 v 6 z ia 2 a c + a + - 
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0. 



(20) 



dW Y AA 2 a (3 Y AA 2 a c 



8(3 



V6 



+ 24V2iY A 2 A bc + M A f3 + 



Y*a 



(21) 
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F «l^ + l) +i -U7! + To71 + To +M *P = - (22) 



4 Analysis 



The purpose of this paper is to give the input parameters appearing in the superpotential 
(H). Though VEVs listed in (||) are functions of the input parameters we will express 
them in the term of the VEVs since we know the desirable values of the VEVs. 



4.1 First Step 

First we check whether it is possible to break SO( 10) down to G2231 consistently with the 
requirement that the spectrum ([!]) remains massless up to 0(e) ~ 0(M VR /Mu)- 

4.1.1 Multiplets under G2231 

First we show what multiplets exist in the SO(10) multiplets (0). 
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I\/Tn 1 "t" l Tilfi't" nnnPT / . , 
iviu.iijijji*3ij u.nu.t-1 *~^2231 


nnH &v**\( ) I 1 (1 i ppitiX i 1 1 oH l n 

U.11U.C1 uv 1 1U 1, UvJlltcLlllCU-lll 


Mfil 

1 N VI X 




^Z, Z, 1 , Uj 


in 19^ 1 9^ 

XU, XZU, XZU 






[ X, 1, O, Zl "T 11. C 


in 1 9fi T% 

XU, XZU, XZU 






(3 i 1 


zl^ 91 n 
^xo, zxu 






fn i nl 

l 1 j O, 1, U J 


zL^i 91 n 
^xo, zxu 




z 


(113 A\ 4- h p 
^x, x, o, — ^-kj ii. c 


zl^i 91 n 
^xo, zxu 


X 


X 


IX, J. , O, U 1 


zt^ 910 

'XU, zxu 






(9 9 3 9^ -U Vi p 
i Zj z, o, zj n. l 


91 n 

^xo, zxu 


y 


y 


(3 i i fiUiip 


1 9R _l l~9fi 

XZU XZU 






(3 1 3 9^ _l V, r 


1 Of\ j_ 1~9fi 

XZU XZU 






(3 1 —9^ _i_ v. p 

1 O, X, U, Zl T 11. L. 


1 oa _l To^ 

XZU XZU 






(1 3 1 R"\ _L h p 
^ x , o , x, — u j ii. l 


1 9fi _l 1~9fi 
xzu xzu 




/C 


(1 3 3 9^ 4- Vi p 
^x, o, o, — zj n. c 


1 9fi _l TOfi 
xzu xzu 






(1 3 ft 9^ -U Vi p 

yLj O , U , Z J 11. L 


1 9fi _l 1~9fi 
xzu xzu 






(9 9 3 A\ 4- Vi p 
^z, z, o, — ^± j n. c 


1 Of, 1 9fi 
xzu, xzu 




- 

y 


(9 9 R -x. Vi p 
i z , z , o , u J n. 


1 9fi 19R 
xzu, xzu 








210 






(1,3,3,-4) + h.c 


210 




X 


(3,1,8,0) + h.c 


210 






(1,3,8,0) + h.c 


210 






(2,2,1,6) + h.c 


210 






(2,1,3,-1) + h.c 


16 + 16 




y 


(1,2,3,1) + h.c 


16 + 16 




X 


(2,1,1,3) + h.c 


16 + 16 






(1,2,1,-3) + h.c 


16+T6 




z 



(23) 



In this table NG1 means a Nambu-Goldstone (NG) mode associated with the break- 
down of SO(10) to Cj223i- An NG mode associated with the SO(10) breaking down to the 
SM group G231 is contained in a multiplet with x, y and z in the column NG2. Under 
G231, certain components of the multiplets with x (y, z) have same quantum number and 
mix with each other. One of combinations of x (y, z) is massless which is swallowed by a 
gauge boson. 

There are also singlets of G2231 which we denote a, b and a. 
4.1.2 F-flat condition 



In the intermediate region c, (3, 
reduced to 



0. And hence the F-term conditions (|T7|) - ( 22 ) are 



dW 
da 



24iV2Y A 2 A ab + M A a = 0, 



(24) 
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(25) 



^ = 24 z V2Y A , A ab - ^ 2Y -f ab + M A a = 0. (26) 
oa 3 

4.1.3 Tuning of parameters 

From now on as we stated at the top of this section, we express the input parameters in 
the term of the VEVs. 

Using the F-flat conditions (E4T) and (KB), M A and M A are expressed as follows: 



M A = M A (Y A 2 A , a, b, a) = ~^^ Y ^ a \ (27) 



M A = M A (Y A 2 A , y AA 2, a, 6, a) = . (28) 

6a 

There is an additional constraint which is obtained by eliminating M A from equations 
and (ESfi: 



- 24 V2 % Y A 2 A a 2 a + AA A -^- + 24 V2 i Y A 2 A ab 2 = 0. (29) 

3 v 2 18 v 2 

We can interpret that this constraint with (P7| ) and (^) is equivalent with that de- 
terminant of the mass matrix for (1,1,3,-4) (= M(l,l,3, — 4) an explicit form is given 
at appendix C) vanishes because (1,1,3,-4) is an NG mode and hence when we substi- 
tute VEVs into the mass matrix for it there must be one massless mode which mean the 
determinant vanishes. 



detM(l,l,3,-4) 

V A M A b Y AA 2 M A b 
M A M A + A 1 A ^ ,_ A 

(30) 



18 3^ 
^ „2 , 1R , V v _ Yl A2 a 2 Y A Y AA2 b 2 



+ 1152Y^2 A a z + 16iY AA 2Y A 2 A aa 



18 108 
= 0. 

Now we required that one (1,1,8,0) mode be massless and therefore determinant of the 
mass matrix for it (= M(l, 1, 8, 0)) should vanish. 
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detM(l, 1,8,0) 

„, Y A M A b Y AA 2M A b 

= M A m a 7= — I 7= — 

18 3^2 

+ ll52Yl 2A a 2 + 16iY AA2 Y A 2 A aa- AA 



(31) 



0. 



108 



Combining fl30|) and (^T|) with substituting (|27|) and (gg), we find 



and 



; i v v 2 , Y A Y AA 2 a a 2 
- y A *a 2 a « H 7^; 1- 16 « Yaa 2 *a 2 a a = 



27 



(((30) - © W& 2 ) 



2304 Y^a 3 + 32 % Y AA 2 Y A 2 A a a 



Y AA 2 a a 1 



Y A Y AA 2 a b 2 



9 54 
-2304 Y^ A ab 2 - 32 i Y AA 2 Y A 2 A ab 2 = 

(((30) + ©)* a). 



(32) 



(33) 



Solving a simultaneous equation (^) and (|33"D we get forms of Y A and Yaa 2 as a 
function of Y A 2 A , a,b,a. Then by substituting these expressions into ( B7p and (B^) we 
find the following three sets of solutions for M A , M4 ,Y A and Yaa 2 as a function of 
Y A 2 Al a,b,a: 




-2A^/2iY A 2 A ab 
a 

2A^2iY A 2 A ab 
-864 iY A 
IUiyI 



2 A a 



(34) 



A 2 A a 



a 
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solution 2: 



solution 3:: 



-2&V2iY A 2 A ab 



A 2 A 



y/2a a 
-432 i Y A 2 A a (-3 



-3b 2 - vV - 10a 2 6 2 + 96 4 
a 2 + 3b 2 - Va 4 - 10 a 2 6 2 + 9 6 4 



(-a 3 + 3a6 2 
-36i Ya 2 a 



a a/o 4 "" 



10 a 2 6 2 + 9 b 4 



aa 



3a 2 + 3b 2 - Va 4 - 10a 2 6 2 + 96 4 



-24 v / 2«>A2 J 4a& 



-24iF A 2 A 6 
\/2a« 

-432 % Y A 2 A a (-3 a 2 + 3 b 2 + Va 4 - 10 a 2 o 2 + 9 6 4 



(-a 2 + 3o 2 + Va 4 - 10a 2 6 2 + 96 4 



-36i Y A 2 A 



+ 3ab 2 + a Va 4 - 10 a 2 o 2 + 9 o 4 



aa 



-3 a 2 + 3 b 2 + Va 4 - 10a 2 6 2 + 96 4 



(35) 



(36) 



In other words, once M A , M A , Y A and Y AA 2 are set to be one of these solutions, the VEVs 
of a, b and a can be chosen at our will and one of (1,1,8,0) mode becomes massless. 

Because we require also that one (1,3,1,0) mode be massless, determinant of the mass 
matrix for it (= M(l, 3, 1, 0)) must be zero. 



detM(l,3,l,0) 
Y A Y AA 2 a 2 1r .,„ ^ Y AA 2® 2 Y A Y AA 2 a b 



— 16 % Y AA 2 Y A 2 A a a 



36 6 lgv ^ 

+ 16 V3iY AA 2 Y A 2 A a b + 1152 Y^ A b 2 + Y ^ M ^ a ( 37 ) 

6 v6 

r- Y A M A b Y AA 2 M A a 

+ 16 Vg i Y A 2 A M A a H — — — AA * +M A M A 

9 v2 v 6 

= 0. 

Using (R7j) and ([54])-(P5]), we obtain a following equation which determine a relation 
between a and b corresponding to a set of above solutions respectively: 
solution 1:: 

a 2 (-3 a 2 + 7 v 7 ^ a b - 6 o 2 ) = 0. 

solution 2:: 

-15 a 6 + 62 V3 a 5 6 + 237 a 4 o 2 - 280 V3 a 3 o 3 - 249 a 2 6 4 + 234 V3 a 6 5 + 27 6 6 
= (33 a 4 - 50 V3 a 3 6 - 78 a 2 b 2 + 78 V3 a b 3 + 9 o 4 ) Va 4 - 10 a 2 6 2 + 9 6 4 . 
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solution 3: 



15 a 6 - 62 VS a 5 6- 237 a 4 6 2 + 280v / 3a 3 6 3 + 249 a 2 6 4 - 234^3 a 6 5 - 276 6 
= (33a 4 - 50 V3a 3 b -78a 2 6 2 + 78 VSab 3 + 9 6 4 ) Va 4 - 10a 2 b 2 + 96 4 . 



Numerically a and b must satisfy the following relation respectively: 
solution 1" 

a - { 



(38) 



solution 2: 



solution 3: 



-0.987293 6 

-0.120361 - 0.724007 z) b 
-0.120361 + 0.724007i) b 
5.112386 



-3.134166 
-0.0643986 6 
1.10047- 0.0616122 i) 6 
1.10047 + 0.0616122 i) b 



(39) 



(40) 



The solution 1 is the exact solution and the others are exact up to O(e). 

In other words, if a and 6 satisfy these relations, one (1,3,1,0) mode becomes massless. 

Other requirements that two (2,2,1,0) modes, one (1,3,1,-6) + h.c mode, one (2,1,3,1) 
+ h.c mode and one (2,1,1,-3) + h.c mode be massless are easily satisfied by tuning 
parameters such as M$, M H , Y H q, A , Y H ^ A and so on. 

To make (1,3,1,-6) + h.c mode massless, from the mass term for it (see appendix C) 



( >/6Y& A a Y^ A a 1$ A 6 N 



(41) 



To make two (2,2,1,0) mode massless we tune parameters M H , M^,Y H ^ A and Y H ^ A 
so that the eigenvalue equation for the mass matrix of (2,2,1,0) 




10 15 J 

+ M9 \ + Y H<5AA*Y H *Ab 2 ^ = Q (42) 

has two solutions (exactly these two solutions may have at most O(e) solution)^. The 
way of getting two zero eigenvalues is to tune the zeroth and first terms of A zero. More 
exactly the zeroth term must be at most 0(e 2 ) and the first term must be at most O(e). 

3 Implicitly it is assumed that the mass matrix for (2,2,1,0) is hermite, that is, all parameters appearing 
in the mass matrix are real 
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To satisfy these constraint 



L]) and the first equation of fl4~3|) lead 



(43) 



\/3g 



<S>A 



6^ 



-K, 



(44) 



up to 0(e). 

Finally to make one (2,1,3,-1) + h.c mode and one (2,1,1,3) + h.c mode massless, foe 
example, we can switch only couplings with subscript 4 on and tune 



Y, 



i/A 



16^ 



iY,p A a/b, 



(45) 



4.1.4 check mass matrices 



%Y^ A a 



7 a 



(46) 



Now we know the necessary condition for the parameters realizing the spectrum ([I]). Then 
we check all the mass matrices to examine whether these parameters really produce the 
spectrum ([!]). 
solution 1" 

The solution 1 does not produce the spectrum (P, because by substituting the solu- 
tion 1 (0) into the mass matrix of (2,2,6,2) multiplet, this multiplet is calculated to be 
massless. 

solution 2:: 

First to see whether the solution 2 , (^) with a relation between a and h Q3Tf), is 
usable, we substitute (|39"D into (j35|). 



M A 
M A 



' 24.3089 «V25a2 a a 
19.1441 iV2Y A 2 A b 2 /a 

" (32.2574 + 5.36258 i) ^2Y& A a 
(-4.78842 + 0.510831 i) V2Y A 2 A b 2 /a 
(-32.2574 + 5.36258 i) V2Y A * A a 
(4.78842 + 0.510831 i) V2Y A 2 A b 2 /a 
-4.69449 iy/2Y A 2 A a 
103.023 i^Y A 2 A b 2 /a 



(47) 
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-13.6527 iY A 2 A b/a 
- , (37.7632 - 7. 13352 z) Y A 2 A b/a . > 

A 12 ^ (-37.7632 - 7.13352 %) Y A 2 A b/a [ } 

677.159 i Y A 2 A b/a 

-104.016 iY A 2 A a /b 

'-1560.23- 131.862 z)y A 2 A a/6 , , 

A * (1560.23- 131.862 i) Y A 2 A a/b [ > 

— 185.139 i Y A 2 A a/b 

In each of these equations, four expressions correspond to the four relations between 
a and b in fl3"9]) respectively. 

As we required that Yukawa couplings are not too big (see the statement below (J?])) 
only the first expression of the solution 2 is meaningful. This means that only the first 
relation between a and b in ( |39| ) is meaningful. 

By substituting (|35"D with the first equation of ( |39l) it is easy to check that all multiplets 
other than those in (|1]) have their mass of O(Mjj) which spread around M\j up to one 
order of magnitude and multiplets in ([!]) are massless. Therefore this solution can be a 
solution of our scenario. 

solution3:: 

First we substitute ( ]4"0"| ) (relation between a and b) into ( |3"ED to see an explicit form of 
solution 3. 



M A 
M A 



Y AA 2 



Y A 



7.65756 iy/2Y A2A a 
-15.8066 iV2Y A 2 A b 2 /a 
372.679 i^/2Y A 2 A a 
1115.98 iy/2Y A 2 A b 2 /a 
(1.21719 - 21.7407 i) V2Y A 2 A a 
(17.3100-22.7812?) V2Y A 2 A b 2 /a 
(-1.21719- 21.7407 1) V2Y A 2 A a 
(-17.3100- 22.7812 i) V2Y A 2 A b 2 / 

—273.079 i Y A 2 A b/a 

3343.29 i Y A 2 A b/a 

(56.3660 + 10.8904 i) Y A 2 A b/a 

(-56.3660 + 10.8904 i) Y A 2 A b/a 

793.766 i Y A 2 A a/b 
6698.93 iY A 2 A a/b 
(241.144- 102.803 i) Y A 2 A a/b 
(-241.144- 102.803 i)Y A 2 A a/b 



(50) 



(51) 



(52) 



In each of these equations, four expressions correspond to the four relations between 
a and b in (|40| ) respectively. 

By the same way as we picked only the first expression up from four cases in solution 
2, the last two relations between a and b in (HO) are meaningful. 
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By substituting fl36|) with the third or fourth equation of (0) it is easy to check that 
all multiplets other than those in (p]) have their mass of 0(Mjj) which spread around Mjj 
up to one order of magnitude and multiplets in ([!]) are massless. Therefore this solution 
can be a solution of our scenario too. 

4.2 Second step 

In this section we find a parameter region which produces our scenario exactly. 
4.2.1 Deviation from the previous solutions 

Because the accuracy of the previous calculation is 0(e), all parameters besides b, a and 
Ya 2 a can deviate from the value which is obtained at the previous section and therefore 
we can expand the deviation in the power of e as follows. 



a = a + 2J a«e* , (53) 

8=1 



M A = M A0 + Y,M Al e\ (54) 
i=i 



M A = M A0 + J2M Ai e\ (55) 

i=l 



Y A = Y m + Y,y^\ (56) 

i=l 



Yaa* = Y AA 2 + ]T Y AAH e\ (57) 
i=i 



/? = EAe\ (58) 



i=i 



c 



£c-^. (59) 



=i 



In these expressions, variables with subscript stand for those which are obtained in 
the previous section. 

Substituting ( |53"D - (|59| ) into the F-flat condition ( p^7|) - (|22|), we get following relations. 
From (0), © and © we get 
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M A i = ai, 

a 

b 3 24y/2iY A 2 A b, b 2 . 

M A1 = — =— Y A1 + ^ ^l + 2-a 1; 60 

9v2a 2 a % 

Vaa^i = (o /6a )Yai H (1 + -ajai- 

a an 



We obtain the relation between /3i and Ci by substituting (|53p - fl5"§D with ( |5Ti| ) into 
(HD and (H) as follows: 

First we note (|T^) and (BTj) can be rewritten 

M(1,3,1,0)(^=-1( 2 ^ A W (61) 

and therefore 



C) 10 V **A 



f)=~M(l,3,l,0)-(^)« (62) 

where M(l,3,l,0) is a mass matrix for (1,3,1,0) and by assumption <f), <fi ~ 0(e). 
Let us decompose the inverse of M(l,3,l,0) 

M(l, 3, 1, 0)" 1 = det (M(l, 3, 1, O))- 1 (A + 0(e)) . (63) 

Since by assumption there is one massless mode in (1,3,1,0) up to O(e), det(M(l,3,l,0)) 
~ O(e) and the first row in A is parallel to the second row in A, that is 

— = — (64) 

a 2l 0-22 

where A = (a^). 



Then up to the leading order of e 



namely, ClS db exact relation 



= a -^c (65) 



A = — ci (66) 



is obtained. 

To see this explicitly, we follow the above calculation in the case of the first relation 
of solution 2. 
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det(M(l, 3, 1, 0)) = (-26423.4^6 a x + 161727iy ^i^ ) £ + Q(e2) 



as we expected the determinant is 0(e). 
A is calculated to be 



A = 

Apparently A satisfies 
Then 



is obtained. 



72.3850 i Y A 2 A ot, -39.5U8iY A 2 A b 
-39.5148 iY A 2 A b, 21.5710 i Y A 2 A b 2 /a 



OL 

fa = -1.83185-ci (67) 
b 



4.2.2 Determination of input parameters of the theory 

Though we can determine the parameters in the power of e order by order, instead of 
doing so we will give the parameters of the theory in a term of the VEVs because the 
purpose of the paper is to find a parameter region for the theory, M's and V's, which 
leads to the spectrum (^). As we will see, by the VEVs a, b, c, a and (3 we can express the 
input parameters of the theory. 

To do this, first we see the F-flat conditions (|i"7f) - (^T|). These equation can be 
rewritten 



where 



C 



( a, 
b, 
c, 
0, 

Vo, 



0. 
0. 
0. 

a, 

ft 



C 



( M A \ 
M A 
Y A 
Y A a 2 

\Y A 2 A J 



? & 2 + 



lSv 7 ^ ' 18^ ' 

^ac+^bc, 



9V2 



0. 

0. 



/i/(ioV6)y*A\ 
i/(ioV2)y* A 
i/ioy* A 

1/5Y<m 



2 V / 6 P ' 

"ST!" 2 ' 



24y/2iab \ 
24V2iaa + 2Ay/2i(3c 
16 VQiac + 24^ibf3 
24V2iab + 8V6ic 2 
2iV2ibc ) 



(68) 



(69) 



As we know from the previous argument that b, c and a can be chosen freely and a 
and (3 are given by 



a + aie, 
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(70) 



where ao is given by the first equation of (|39"1) or one of the last two equation of (|40D and 
/?i is given by fl66|) . Note that higher orders in d53|) and (|58"D can be absorbed into a x and 
/3 2 respectively. 

Then the input parameters are reduced to 



detC 



/ M A \ 


( 


M A 




Y A 


= -c- 1 






\Ya*a) 


V 


case of solution 2, 


(det C) 





i/(ioV2)y$ A 

1/IOF.ka 

\/6/ioy^ 
i/5r*A 



(71) 



3.76350 i a 2 6 4 /3 2 c a - 2.25347kt 3 6 2 aiC 2 ) e 3 + 0(e 





( > 


0. 


-2.68018 i a 3 b 3 c u 


0. 


-1.08826 i a 4 b 2 d \ 




o, 


0. 


-2.11074 % a b 5 c u 


0. 


-0.857040 % a 2 b A c x 




o, 


0. 


8.10927ia 3 b 2 c u 


0. 


3.29268 z a 4 6 ci 




o, 


0. 


1.06439 i a 6 4 Ci, 


0. 


0.432184 i a 2 6 3 ci 




Vo, 


o. 


-0.0779620 a 2 b 3 c h 


0. 


-0.0316556 a 3 b 2 cj 


+ 0(6) 











From this equation it is easy to see that all parameters are of order e° and they satisfy 
the first solution of the solution 2. 

Finally from (|2~2"D M$ is determined: 

= - y « br + f J " y - (iwi + wf + To) ■ (72) 



4.2.3 check mass matrices 

The multiplets in ([I]) besides one (2,2,1,0) must decouple at M^that is, they must acquire 
mass of 0(M VR ). 

From now on we check whether they have mass of 0(M UR ). 

First we note one (2,1,3,-1) + h.c and (2,1,1,3) + h.c can have mass of 0(M UR ) by the 
following two reasons: (1) Parameters Y^a and may deviate from the value given by 
( f45|) and fl4"6| ) respectively Q. (2) There exist couplings with c and 3. 

4 Though (2,1,3,-1) + h.c has a same quantum number under the SM group as an NG mode associated 
with the breakdown of SO(10) the SM group (see table (|23|)), it does not mix with others because the 
VEV of ijj = and therefore this NG mode does not consist of it. (2,1,1,3) + h.c has a same quantum 
number as that of (2,2,1,0) under the SM group but by the same reason they do not mix with (2,2,1,0). 
See the superpotential (||) - (||). 
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Then we see the mass matrix for (2,2,1,0). Under SM it has a quantum number 
(2,1, ±1/2). (2,2,1,6) + h.c also includes the same component. Then the mass matrix is 



where 



M(2, l,±l/2) 



fM A , 


X, 


y- 





\ 


x', 


M H , 


u, 


V 




o, 


u, 


o, 


,w- 


z 


V y', 




w + z, 





J 



M A 
x 
x' 

y 
y' 

u 

V 

w 

z 



M(2, 2, 1, 6) + ±Y A c + 24iY AA2 p 
-^W>~0(e 3 / 2 ) 
-^Y mA <j> ~ 0(6 3 / 2 ) 



-iy*A0 ~ 0(e) 



— tkY$, A 



40 

—k*Y, 



<p ~ 0(e) 



y$A c _|_ 



30 



10 



O(e) 



(73) 



(74) 



M(2,2,l,6) is given in the appendix C. Orders of x, y, .. are followed from fl43|) 

Because one (2,1, ±1/2) multiplet remains massless after G2231 breaks down to the SM 
group 



det (M(2, 1, ±1/2)) = {M A (z 2 - w 2 ) +yy'(w - z)}M H + 2M A uvw. 
and hence M H is determined as follows: 



0. 



(75) 



2uvw ^ . . 
M H = — -±0(e). 



■ur 



(76) 



In this case the higher order terms must be included to have a pair of light Higgs doublets. 

Next let us consider (1,1,8,0). This multiplet becomes (1,8,0) under the SM group 
and therefore it mixes with T 3 # = component of (1,3,8,0) under the SM. Then the mass 
matrix for (1,8,0) is represented as 3 x 3 matrix. 



Af (1,8,0) 



M(l, 1,8,0) 


mixing 


mixing 


M(l,3,8,0) 



(77) 



After G2231 breaks down to the SM group, there is a correction of 0(M;ye ~ M VR ) to 
the mass matrices M(l, 1,8,0) and M(l,3,8,0) because parameters appearing in them are 
different by O(e) from those calculated at the previous section. It is directly calculated 
using (fflf) (or equivalently ( |53D - ( [j7f) and (|60|) ) that one of the eigenvalues of M(l, 1,8,0) 
is of O(Mjj) which has already suggested at the previous section and the other is 0(M UR ). 
As M(l,3,8,0) is O(Mu), even though there is a correction of 0(M UR ), M(l,3,8,0) is still 
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O(Mu). Contributions of c and (3 to the mass matrix (|77|) appear at mixing terms between 
(1,1,8,0) and (1,3,8,0)0 and they are of 0{M VR ). Then M(l,8,0) takes the following form 

0(M ue ) 
0{Mue) 0(M ue ) I . (78) 




0(M ue ) O(Mu) 

Apparently two eigenvalues are of O(Mjj) and the other is of 0(M UR ). This fact 
suggests that the lightest element of (1,1,8,0) under G2231 decouples at the scale M VR . 

Finally we check the mass of (1,3,1,0) and (1,3,1,-6) + h.c. Under the SM (1,3,1,0) 
is decomposed into one neutral singlet and a pair of charged singlet with hypercharge 

Y = ±1. (1,3,1, ,-6) + h.c becomes two neutral singlets, a pair of Y = ±1 and a pair of 

Y = ±2 singlets. Then Y = ±1 component of them will mix with each other. 

Mass for Y = ±2 component takes the following form 

Y - " I ) + r « (wJS + - + M * = -l Y ^ - (79) 



where (f72|) is used. 

From this equation obviously the Y = ±2 component has a mass of 0(M UR ). 
Mass matrix of Y — ±1 component is 

/ -X±Ap + M A , _2^ + 2 4z v / 21a2 A 6, 

-L^L + 2 Aiy/2Y A 2 A b, f^ + WiVSY A * A a+^l + M A , 

y<s>A <t> Yj'A <t> y<s>A p _ g 

5 ' 10 ' 5 10 ' , , 

(80) 

Since it is an NG mode associated with the breakdown of G2231 to G231 there is one 
massless mode. It is easy to see that this matrix has eigenvalue because 1st row x (3/4> + 
2nd row x c/0 + 3rd row = using the F-flat conditions (|19"D and (0). It is also explicitly 
calculated that one eigenvalue is of O(Mu) and the other is of 0(M UR ). 



5 summary 

As we saw, by constructing the input parameters for the theory using (|7l|), (|72"D, (|74|) and 
( |76|) from the desired values of VEVs a, b, c, a, /3, and (ft which satisfy (|T0"D and (|70|), we 
can have particles ([^) in the intermediate region. They decouple from the spectrum at 
M VR except a pair of what we call Higgs doublets. 

It means that it is possible to construct a SUSY SO(10) GUT with an intermediate 
scale consistent with the gauge unification. It suggests also that the right-handed neu- 
trinos acquire mass through a renormalizable coupling, and it can be understood as a 
reflection of the breakdown of G2231 to G 2 3i 

5 There is no contribution of c and f3 to M(l, 1,8,0) and M(l,3,8,0). The reason is as follows. Under 
G2231 c and P are contained in (1,3,1,0). Because (1,3,1,0)(1, 1, 8, 0) 2 contains no singlet, c and (5 do 
not couple to (1,1,8, 0) 2 . Though (1,3,1,0)(1, 3, 8, 0) 2 can appear, as there is no three point coupling of 
T3R = component of SU(2) triplet, c and (3 do not couple to T^r = component of (1,3,8,0) 
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There are many variations for a SUSY SO(10) GUT with an intermediate scale because 
there are many candidates for the particle content which exist in the intermediate region 
and we have many variations for content of SO(10) multiplets which contain one of the 
candidates. 

For example, we can replace (2,2,1,0) to (2,1,1,3) + h.c in the spectrum (|l|) and vise 
versa, because their contribution to the running of the gauge coupling relevant to G231 is 
as same as that of each other. 

When we remove one (2,2,1,0) from the spectrum ([]]) and add one (2,1,1,3) + h.c to it, 
by adding a pair of SO(10) multiplets 16 + 16 which contains (2,1,1,3) + h.c under G2231 
we can have such a spectrum at the intermediate region. At that time while we have to 
tune couplings relevant to SO(10) multiplets 16 + 16, we can release the constraint ( f43|) 
(or equivallently (|74f )). 

Of course, there is a quite different type of content for the candidates. Using them we 
can construct quite a different SO(10) GUT with an intermediate scale. 

Though the gauge unification by the MSSM is a very attractive idea,to take into 
account a right handed-neutrino mass we should consider a possibility of a GUT with an 
intermediate symmetry. 
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A The reason why we need a multiplet (1,3,1,0) 

Here we show the reason why we need a multiplet (1,3,1,0) in the intermediate region. 

First we note that we required at least there is a pair of multiplet (1,3,1,-6) + h.c 
(= <3> + $) in the intermediate region || and hence at this region in the superpotential 
effectively there must be a term 

W = M $ $$. (81) 

Because we consider an SO (10) GUT the mass parameter M$ is ,in general thought 
to be of 0(Mu). 

In this case it is, however, impossible that $ acquires a VEV. Of course if we tune the 
parameter M$ be 0, as there is a flat direction in D-term, $ can acquire a VEV, but in 
this case there are two problems: 

(1) there is no way to determine a magnitude of the VEV of $. 

(2) hypercharge Y = ±2 component of $ cannot have massQ. 

6 Note that only an NG mode can get a mass through D-term. In general, such a component corresponds 
to a massive gaugino. 



22 



Then we have to add other multiplets. The easiest way to solve the problem (1) is to 
add a singlet (= S) []. If there is a singlet the superpotential will have a form 

W = M $ $$ + Y 9S S^ + \m s S 2 + ^Y S S 3 (82) 
and F-flat conditions are (< $ >= (f),< S >= s) 

dW - 

— - = (M* + r„ s s)0 = 0, (83) 



dW - 1 

— - = F $s # + M s s + -Y s s 2 . (84) 
as 2 

Then VEVs are determined to 



Ki 



*5 



(85) 



Y$S 1 

Though as we mention below (^l|) M's are thought to be of O(Mjj), we can give a 
VEV of 0(M UR ) to $ if coupling constants are fine-tuned while s is of 0(Mu). 

Unfortunately even after we add a singlet, the problem (2) is not solved because the 
mass for Y=±2 component is 

M $ + Y s s = (87) 

according to the F-flat condition (^3|). The reason why it is still massless is that no 
multiplet couples to $ which acquires a VEV of 0(M UR ) and distinguishes the component 
of a SU(2) R triplet and hence all component of $ is still degenerate after SU(2) R breaking. 

This means that to make Y — i2 component decouple from the spectrum after SU(2)^, 
breaking we have to make a multiplet couple to $ which will get a VEV of 0(M VR ) and 
distinguishes the component of a SU(2) R triplet, that is, a non-singlet. It is easy to find 
what non-singlet can couple to $$. From we have three representation: 



(1,1,1,0) 

(1,3,1,0) (88) 
(1,5,1,0) 

As SU(2) fl non-singlets are the latter two and (1,5,1,0) is not contained in a relatively 
smaller representation of SO(10), we have to use (1,3,1,0). Since T 3R = component of a 
triplet is an SM singlet it can get a VEV. 

7 Because we consider an SO(10) GUT, there are several singlets though naturally their masses are of 
0(Mu). 
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Since (1,3,1,0) is not a singlet under G2231, its VEV is at most of 0(M UR ), while 
because (1,3,1,0) gives a mass of 0(M UR ) to Y=±2 component of $, even if there are 
many (1,3,1,0), one of their VEV must be of 0(M UR ). This implies that at least one of 
(1,3,1,0) must have a mass of 0(M UR ). In the following we will see it explicitly. 

First when there are also (1,3,1,0) multiplets (= B{) the superpotential takes a follow- 
ing form. 

W = M$$$ + YisSQQ + Y, Y i R i^ 

i 

+ Im s S 2 + ^Y s S 3 (89) 

and F-flat conditions are (< Bi >= fa) 

dW 



(M» + r* s s + x;y i /3 < )$ = o, (90) 



° ] 1 W0 + M s s + l -Y s s 2 + £ Y Sij fafa = 0, (91) 



— - = YM + £(My + V ijS )A = 0. (92) 

Note that there is no three point coupling of T 3 = component of SU(2) triplet and 
hence there is no affect of Y^. 
From (|92|) fa is calculated to 

fa = -{M- 1 )^^, (93) 
M i:j = (Mij + Yjs). 

By assumption (f> ~ 0(M VR ) and as we mentioned one of fa also must be of 0(M VR ). 
These facts imply that in the above equation M must have at least one eigenvalue of 
0(M UR ). Because M is a mass matrix for (1,3,1,0) (see (|S9"D), it means that at least one 
of (1,3,1,0) must be massless at the GUT scale. 

In this case mass for Y= ±2 is calculated 

(M* + Y* s s - J2 = "2 E °*Pi ~ °( M ^) ( 94 ) 

where (|90f) is used. Apparently this component decouples at M UR , namely, the problem 
(2) is solved. 
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B Construction of Representations 



In this section we briefly review how we construct representations of subgroups contained 
in SO (10) representations and give the rule for calculating CG coefficient appearing in 
three point couplings. However, we do not mention about an SO(10) spinor 16 because 
it is impossible to understand the meaning of the indices for a spinor in the same way 
of understanding that for an SO(10) vector 10 and essentially we do not need to handle 
them directly in this paper. To see how to handle an SO(10) spinor, see ref. ||12|| . When 
calculating CG coefficient relevant to a spinor the gamma matrices for SO (10) constructed 
explicitly in the reference are used. 

B.l Meanings of Subscripts 

For SO(10) the fundamental representation is a 10 dimensional real vector 



H = (Hi), * = 1, ...,10. 

It means when we construct a fundamental representation for SO(10) we can use a 
following basis for it: 



H = hie i: 



(95) 



where 



Hi 6^ H 'j 6j 





1 





}i th component. 



(96) 



Here after in this appendix, repeated subscripts are assumed to be contracted. 
In this case index % means nothing but SO(10) vector. 

For our convenience we can attach an additional meaning to it. SO(10) includes SU(5) 
<8> U(l) and SO(6) <g> SO(4) ~ SU(4) ® SU(2) ® SU(2). Under them the fundamental 
representation 10 is decomposed into [O] 



5(2) + 5(-2) under SU(5) <8> U(l) 

10 = { (6, 1) + (1,4) under SO(6) g> SO(4) 
(6, 1,1) + (2, 2, 1) under SU(4) ® SU(2) 



SU(2) 



^Exactly in a mathematical term what fundamental representation means is identity representation. 
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Then we can add a meaning of, for example, SO (6) vector to indices 1 to 6 and SO (4) 
vector to 7 to 10[]. Here after stands for 10. In other words SO(6), an SO(10) subgroup, 
acts on the indices 1 to 6 and SO(4) acts on 7 to 0. 

We can add more meaning to indices of an SO (10) vector by giving a meaning 5(2) 
representation under SU(5) <g> U(l) to (1 + 2i, 3 + 4i, 5 + 6i, 7 + 8i, 9 + 0i) and its complex 
conjugate to (1 — 2i, 3 — 4i, 5 — 6i, 7 — 8i, 9 — Oi). 

What 1 + 2i means is as follows. When we construct a vector representation we can 
use a basis E a+bi and its complex conjugate E a -u = E a+ bi where b = a + 1 and a is an 
odd number other than which is introduced at the top of this section. 



E 



a+bi 



V2 



( \ 
1 

i 

V o J 



}ath 1 i 

}bth = 7! ea + 7! efc 



(97) 



where -j= is a normalization factor to achieve E\ +bi E a+bi = 1. 
Then 



H — hiei — h a+ biE a+ bi + h a -uE, 



a—bi 



where 



E l+u H 



V2 



(h a - h b i) 



(91 



h a+ bi is a component of a SU(5) vector and its U(l) charge is 2. As it is easily seen the 
component for an SO(10) vector depends on a basis. 

Because both SU(5) and SO(6) ~ SU(4) contain SU(3) C we can add the meaning of 
SU(3) 3 and 3 to the SO(6) vector indices 1 to 6: (1 + 2i, 3 + 5 + Qi) is an SU(3) vector 
3. By the same way we can add the meaning of SU(2) 2 and 2 to the SO(4) ~ SU(2) ® 
SU(2) vector indices 7 to 0: (7 + 8i, 9 + Oi) is an SU(2) vector 2. 

As we will see later a higher representation is represented as a tensor. By this con- 
struction when we consider what representations a higher representation contains under, 
for example, SO(10) subgroup SU(4), it is sufficient to deal with indices 1 to 6. When 
considering SU(5) subgroup we can deal with combinations of SO(10) subscripts 1 + 2i 
and so on. 



B.2 SO(10) Representations and Representations of subgroups 
contained in SO(10) Representations 

The representations 45, 126 + 126 and 210 are formulated from the fundamental repre- 
sentation as antisymmetric tensors of 2nd, 5th and 4th rank respectively. By the charac- 
teristic of SO(10) 5th rank antisymmetric tensor is decomposed into two parts, 126 and 

9 In the papers |j| ^] they give a meaning of SO(6) vector to indices 5 to 10 and that of SO(4) to 1 to 

4 
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126. Using 10th rank antisymmetric e tensor (= e a b c deijkim) it is decomposed into two 

"^^■abcdeijklm^ ijklm ~^~^abcdei 



eigenstates 12 



£ abcdeij klm® ij klm ^ abcde- 



(99) 



What has a plus eigenvalue is defined to be 126 and the other is to be 126. 

In the same way as an SO(10) vector 10 we can express these representations using a 
component and a basis. To express 45 (= A) we can take a basis ey as follows: 

A = dijtij (100) 

where 

aij = tvAeij, e y = ((ey)ab) = -y=(6 ai 5 bj - 5 aj 5 bi ). (101) 

ciij corresponds to a component of 45 representation. In our notation subscripts i,j for a 
component and a basis satisfy that i > j. 

In a similar manner 126 + 126 (=$ + $) is written as 

$ (or $) = (pijkimeijkim (102) 

where eijkim is an antisymmetric tensor and only when a combination of indices coincide 
with subscripts {ijklm} it has a value 1/V5! or — l/v5! . The sign is defined to make 
Zijkim be antisymmetric. Here {ijklm} satisfies i>j>k>l>m. Exactly, for e^m 
to be a basis of 126 (or 126) there is another constraint for it as we explained at (p9|), 
though we do not touch the detail here. Then a component of 126 is given by 



'ijklm 



&abcde( &ijklm) abcde- (103) 



is a necessary normalization factor to express a 126 representation by (|102|) and (|103|) 
similar to ^= in (|97D . 

In the case of 210 a basis for it becomes 4th rank antisymmetric tensor and its nor- 
malization is l/v^4L Besides it 210 (= A) is represented in the same way: 

A = dijkieijki 

where 

&ijkl ^abcdi&ijkl) abed 

and i > j > k > I. 

To construct a representation under subgroups we use a linear combination of these 
basis in the same way that when we extract a 5(2) of the subgroup SU(5) ® U(l) from 
an SO (10) vector we use a basis E a+bi . 
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For example let us consider G231 singlets contained in 126 and 126. They are SU(5) 
singlets. Then it is sufficient to deal with SU(5) subscripts 1 + 2? and so on. By the 
quintality of SU(5) the form of the basis of SU(5) singlets in 126 and 126 are determined 

to be ei_2i, 3 -4i,5-6i,7-8i,9-0i,ei+2i,3+4i,5+6i,7+8i,9+0i- They are understood in the same way 
as E 1+2l , 



ei-2i,3-4i,5-6i,7-8i,9-0i 



(<3l3579 — ^23579 + •••) 



where h= is an extra normalization factor to achieve 



( e l-2i,3-4j,5-6i,7-8j,9-0j)a6cde( e l-2j,3-4i,5-6j,7-8i,9-0j)a6cde 



similar to -4= m 



1 

V2 



It is easily seen that the former is a basis of 126 and the latter is that of 126 by making 
Sabcdeijkim acting on them or by counting U(l) charge fTTH. All other representation of 
subgroups contained in SO(10) representations are constructed in a similar way. 



B.3 CG coefficient 

Using 10, 45, 126, 126 and 210 we have following SO(10) singlets |TJ 

#$A, H$A, A 3 ,$A$, <M$, A 2 A, AA 2 
We can get singlets by contracting all indices of tensors: 



H$A 


— H a ^ ' abode A^cde 


Hl>A 


= Ha&abcde Abode 


A 3 


= A abcd A cde fA e f ab 


$A$ 


= ^abijkA a f)cd^cdijk 


<M$ 


= QaijklAab&bijkl 


A 2 A 


= ^-ab^-cd^abcd 


AA 2 


= & abcdefghijA a \)A c( ief A 



ghij 



In a term of components of the representations 



H$A 
H$A 
A 3 
$A$ 



ha<Pabcde $bcde , 

^^y^^abcd^cdef^efabi 
1 - 



abijk Vabcd4 ) cdij k > 
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<t>A<t> = ^aijklAab^Ujkl, 



5\/2 

1 

abcdi 



A 2 A = -—^a ab a cd 5, 

AA 2 = 24V2ia ab 5 cdef 5 ghij . 

where repeated subscripts are not summed and in the last equation abode fghij are dif- 
ferent from each other. 

Then we rewrite the superpotential (|5|) in a term of components, for example, 



Y&A 3 — - ^ 5 a bcd5cdefS e fab 



and so on. Therefore for components that as an expansion parameter for the perturbation 
Yukawa coupling = 1 means = 6v^6 and so on. 

Of course, since a component of an irreducible representation is a linear combination 
of these components, CG coefficient for an irreducible representation is different from, for 
example, ^? in the case of A 3 . 

For example let us calculate a CG coefficient for the singlet (3 contained in 45 and a 
contained in 210 ( see the table (|3])). They are contained in the form A 78+90 = /3e 78+9 o 
and A 7890 = ae 7890 respectively. Then 

A a bA c( iA a b c d = P 2 a(ejs+9o) ab(^78+9o) cd(^789o) abed 

In the second line | comes from an element of e7§ + go and comes from an element of 
e7 890 . 2! comes from a summation between {ab} and {cd}. {ab} and {cd} are {78} or 
{90}. The last factor 2 comes from an exchange of {78} and {90}. 



C Mass matrices under G2231 and their eigenvalue 
equations 

Under G2231 the multiplets of our model have mass terms as follows. They are listed 
following the order of the list (|23|). Full mass matrices are given with contributions 
from c,/3,(f> and after G2231 breaks down to G23i- But these contributions are of order 
M UR ~ Mjje and hence if the mass eigenvalue is of O(Mjj), they are negligible and we do 
not need to consider them. 
(2,2,1,0) multiplet; 
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M(2,2,l,0) = 
'1,1,3,-2) + h.c multiplet; 



Yh<s>a b 



V 



M(l,l,3,2) = 



( M H , 

Yh<s>a (y/Sa—b) 



30 



VTo ' 



V10 ' 



V10 

y% A b + M 



VlO ' 15%/2 



°> 1573 



Yhza (vgo-b) V g¥A (v^a+fe) \ 
v/30 ' v/30 









(3,1,1,0) + h.c multiplet; 



y AA2 a 



M(3, 1,1,0) = 
(1,3,1,0) multiplet 



■ ^p-24iV2>W> 



y AA2 a I JIT 

M(l,3,l,0)- I ^ + M * vTT 



^A^ + 24iv / 2lA^& 



(1,1,3,-4) multiplet; 

M(l,l,3,-4) : 
(1,1,8,0) multiplet; 



Y Ab , jit 
18 + MA 



3V2 



2AiV2Y A 2 A a-^^, 



M(l, 1,8,0) = 
(2,2,3,2) + h.c multiplet; 



3^ 



24iy/2Y A 2 A a-?fff' 



24^^^-^^, 
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/ M A , 



M(2,2,3,2)= Sv^^A 



3 ' 



y AA2 ° 



8y / 6iF A 2 j4 6, 

M A , 16i"N/3y A 2AQ! 
16^^ a, ^= + M A 



30 



(3,1,1,6) + h.c multiplet; 



M(3, 1, 1, 6) = — 7= H 7= + M$ 

v ; 10 10 10 V2 



(3,1,3,2) + h.c multiplet; 



13, 1, 3, 2 = ^= - -^-= + + M$ 

5^ 10 V6 30^ 



(3,1,6,-2) + h.c multiplet; 

MC\ 1 fi 9^ - Y * a0l Y * Aa Y ^ b 4. A/f 



(1,3,1,-6) + h.c multiplet; 



M(l, 3, 1, -6) = — 7= H 1= + M$ 

V 7 10 10^ 10\/2 



(1,3,3,-2) + h.c multiplet; 

M(l,3,3,-2) = ^ + ^ + ^ + M* 
5^ 10 V6 30 V2 



(1,3,6,2) + h.c multiplet; 



M(l, 3, 6, 2) = Y$Aa _ Y$&b 

5^6 10^6 30^ 



(2,2,3,-4) + h.c multiplet; 
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M(2,2,3,-4) 



15 + 30V2 +M *> U 

U ' 15 + 30^ + iW * 



(2,2,8,0) multiplet; 



*(2,2,8,0) = -^ + tf. 



(3,1,3,-4) + h.c multiplet; 



M(3, 1, 3, -4) = -|^= -8^^^+ ^= + M A 



(1,3,3,-4) + h.c multiplet; 



M(l, 3, 3, -4) = ^ + 8 ? >/6y A »A « + ^= + M A 



(3,1,8,0) multiplet; 



M(3, 1, 8, 0) = + 8 i Ve Y A , A a - ^= + M A 



(1,3,8,0) multiplet; 



M(l, 3, 8, 0) = ^ - 8 i V6 Y A2A a - ^= + M A 



(2,2,1,6) + h.c multiplet; 
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M(2, 2,1, 6) = M + Ma 



(2,2,6,-2) + h.c multiplet; 



M(2,2,6,-2) = --^4 + M A 



(2,1,3,-1) + h.c multiplet; 



M(2, 1,3,-1) 



(1,2,3,1) + h.c multiplet; 



18 







-yY^Aa + 2^ 3A ( v / 6a + v^) 
V -^«^ 44 a + 2F v , 4A ( v / 6a + v 7 ^) + My j 



M(l,2,3,l) 



^il^Asa + 2y* A3 (-v / 6a + v^) 
V ^giy*A4a + 2y tA4 (-v / 6a + v 7 ^) + My J 



(2,1,1,3) + h.c multiplet; 



M(2,l,l,3) 
(1,2,1,-3) + h.c multiplet; 



( r" N 

V6iYy A3 a + 2v / 6Y^ A3 (a — V3b) 
V V6iYy A4 a + 2 v / 6y* A4 (a - v 7 ^) + 7 



M(l,2,l,-3) 



/ 

— y/EiY$A2<X 

— y/6iYy A3 a — 2^/6Yy A3 (a + V3b) 
V -VGiYvAta - 2V6Yy Ai (a + V3b) + My J 
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